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' The effects of generalized uncertainty principle (GUP) on the inflationary dynamics and the 

. thermodynamics of the early universe are studied. Using the GUP approach, the tensorial and 

scalar density fluctuations in the inflation era are evaluated and compared with the standard case. 
We find a good agreement with the Wilkinson Microwave Anisotropy Probe data. Assuming that 
<n: a quantum gas of scalar particles is confined within a thin layer near the apparent horizon of 

the Friedmann-Lemaitre-Robertson- Walker universe which satisfies the boundary condition, the 
' number and entropy densities and the free energy arising form the quantum states are calculated 

using the GUP approach. A qualitative estimation for effects of the quantum gravity on all these 
I ^ ■ thermodynamic quantities is introduced. 
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If} . I. INTRODUCTION 

The idea that the uncertainty principle would be affected by the quantum gravity has been suggested couple 
decades ago Q . Should the theories of quantum gravity, such as string theory, doubly special relativity and black 
hole physics be confirmed, our understanding of the basic laws and principles of physics turn to be considerably 
different, especially at very high energies or short distances 04^] • Various examples can be mentioned to support 
this phenomena. In the context of polymer quantization, the commutation relations are given in terms of the 
polymer mass scale Q. Also, the standard commutation relations in the quantum mechanics are conjectured 
to be changed or better to say generalized at the length scales of the order of Planck's length HI. Such 
modifications are supposed to play an essential role in the quantum gravitational corrections at very high 
energy . Accordingly, the standard uncertainty relation of quantum mechanics is replaced b y a, gravitational 
uncertainty relation having a minimal observable length of the order of Planck's length (6l. llXj| - jl3| . 

The existence of a minimal length is one of the most interesting predictions of such new physics. These can be 
seen as the consequences of the string theory, since strings can not interact at distances smaller than their size 
which leads to a generalized uncertainty principle (GUP) 0. Furthermore, the black hole physics suggests that 
the uncertainty relation should be modified near the Planck's energy scale because of measuring the photons 
emitted from the black hole suffers from two major errors. The first one is the error by Heisenberg classical 
analysis and the second one is because the black hole mass varies during the emission process and the radius of 
the horizon changes accordingly 0, 0, [HHItJ . As discussed, these newly-discovered fundamental properties of 
space-time would result in different phenomenological outcomes in other physical branches (lcij . In the first part 
of this present work, we want to investigate the effects of GUP on the inflationary parameters in the standard 
inflation. 

At very short distances, the holographic principle for gravity is assumed to relate the gravitational quantum 
theory to quantum field theory. At this short scale, the entropy of a black hole would be related to the area 
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of the horizon pjl, [2(| ■ The covariant entropy bound in the Friedmann-Lemaitre- Robertson- Walker (FLRW) is 
found to indicate to a holographic nature in terms of temperature and entropy [2l[ . The cosmological boundary 
can be chosen as the cosmological apparent horizon instead of the event horizon of a black hole. In light of 
this, we mention that the statistical (informational) entropy of a black hole can be calculated using the brick 
wall method [22I ]. In order to avoid the divergence near the event horizon, a cutoff parameter would be utilized. 
Since the degrees of freedom would be dominant near horizon, the brick wall method is used to be replaced by a 
thin- layer model making the calculation of entropy possible |23l - [30( . The entropy of the FLRW universe is given 
by time-dependent metric. The GUP approach has been used in calculating the entropy of various black holes 
[3ll - l42j . The effect of GUP on the reheating phase after inflation of the universe has been studied in [43|. The 
present work aims to complete this investigation by studying the effect of GUP in the inflationary era itself. In 
doing this, we start from the number density arising from the quantum states in the early universe. Then, we 
calculate the free energy and entropy density. The idea of calculating thermodynamic quantities from quantum 
nature of physical systems dates back to a about one decade [3j49j , where the entropy arising from mixing of 
the quantum states of degenerate quarks in a very simple hadronic model has been estimated and applied to 
different physical systems. 

Some basic features of the FLRW universe are given in section [TTJ The GUP Approach which will be utilized 
in the present work is elaborated in section Mil The whole treatment is based on the inflation era. Section IIVI 
is devoted to the second topic of the present work, some thermodynamic quantities arising from the quantum 
states of the early universe. The conclusions are listed out in section fVl 



II. THE FLRW UNIVERSE 



In the FLRW universe, the standard (n + l)-dimensional metric reads 

ds 2 = h ab dx a dx b + r 2 dn 2 n _ x , (I) 

where x a = (t,r) and h a b = diag(— 1, a 2 /(I — kr 2 )). d$l 2 l _ 1 is the line element of an n + I-dimensional unit 
sphere. a(t) and k are scale factor and curvature parameter, respectively. Then, the radius of the apparent 
horizon is given by 
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Ra =[ H+ ^) • < 2 > 

It is obvious that the time evolution of the scale factor entirely depends on the background equation of state. 
Seeking for simplicity, we utilize (50j 

a(t) = t 2 / 3k , (3) 

where t is the cosmic time and k = 1 — {be) 2 /(I — c 2 ). The parameters b and c are free and dimensionless. 
Their values can be fixed by cosmological observations. Then, the Hubble parameter and radius of the apparent 
horizon read 

o 1 

(4) 
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(5) 

From the metric given in Eq. (TTJ) and the Einstein in non-viscous background equations, we get 
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H+ ^ = — p+ r (6) 



= -^G(p + P ). (7) 
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Then, the total energy p and temperature T inside the sphere of radius Ra can be evaluated as follows. 
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where n gives the dimension of the universe. From Eq. ^ and ([5]) , it is obvious that the inverse radius of the 
apparent horizon is to be determined by the energy- momentum tensor i.e., matter and cosmological constant 
A. G is the gravitational constant and p is the pressure. Taking into consideration the viscous nature of the 
background geometry makes the treatment of thermodynamics of FLRW considerably complicated [5ll - l59| . For 
completeness, we give the cross section of particle production 
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where Y is the gamma function and M p i is the Planck mass. 



III. TENSORIAL AND SCALAR DENSITY FLUCTUATIONS IN THE INFLATION ERA 



At short distances, the standard commutation relations are conjectured to be changed. In light of this, a new 
model of GUP was proposed [60l - l62j . It predicts a maximum observable momentum and a minimal measurable 
length. Accordingly, [xj,Xj] = [pi,Pj] = (via the Jacobi identity) turn to be produced. 
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is Planck's mass 



where the parameter a = ao/M p c = ao£ p /h and M p c 2 stand for Planck's energy. 

and length, respectively, ao sets on the upper and lower bounds to a. Apparently, Eqs. (fTTj) imply the existence 
of a minimum measurable length and a maximum measurable momentum 
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where Aa; > Ax m i n and Ap < Ap max . Accordingly, for a particle having a distant origin and an energy scale 
comparable to the Planck's one, the momentum would be a subject of a modification [60l462l ] 
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where X{ = xoi and poj satisfy the canonical commutation relations [xoi,poj] = ^ $ij an d simultaneously fulfil 
Eq. (|lip . Here, poi can be interpreted as the momentum at low energies (having the standard representation 
in position space i.e., poi = —ihd/dxoi) and p, as that at high energies. 

As given in [63| and Eq. (fTTj) . the first bound for the dimcnsionless ao is about ~ 10 17 , which would 
approximately gives a ~ 10~ 2 GeV -1 . The other bound of ao which is ~ 10 10 . This lower bound means that 
a ~ 10 -9 GeV -1 . As discussed in [gj], the exact bound on a can be obtained by comparing with observations 
and experiments f65| . It seems that the gamma rays burst would allow us to set an upper value for the 
GUP-charactering parameter a. 

In order to relate this with the inflation era, we define <fr as the scaler field deriving the inflation in the early 
universe. Then, the pressure and energy density respectively read 



P{4>) = l<j> 2 -V(4>), 
p{<t>) = \i> 2 + v{<p), 



(15) 
(16) 



where V(<p) is the inflation potential, which is supposed to be sufficiently flat. The main potential slow-roll 
parameters [rP^ ] are given as 
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« = M *vwr (18) 

where M p i = m v i/"J&K is a four dimensional fundamental scale. It gives the reduced Planck's mass. The 
slow-roll approximations guarantee that the quantities in Eq. (|17[) and (jf 8[) are much smaller than unity. These 
conditions are supposed to ensure an inflationary phase in which the expansion of the universe is accelerating. 
The conformal time is given as 

where a is the scale factor and H = a/a is the Hubble parameter. 

In order to distinguish from the curvature parameter fc, which is widely used in literature, let us denote the 
wave number by j. Here, j is assumed to give the comoving momentum. It seems to be r-dependent and can 
be expressed is terms of the physical momentum P 

3=aP = ~ (20) 

In the GUP approach, the momentum is subject of modification, j — > j(l — aj). Accordingly, the modification 
in the comic scale a reads 

a=^l. (21) 

Then, in the presence of the minimal length cutoff, the scalar spectral index is given by 

d lnp s . , d lnp„ , , 

amj(l-Qj) a In j 

where p s is the amplitude of the scalar density perturbation i.e., the scalar density fluctuations. Due to the 
modified commutators, a change in H is likely expected. This can be realized using slow-roll parameters. In 
the standard case, the spectral index can be expressed in these quantities [66| . 

n s = l + 277-6e, (23) 

where rj and e are given in Eqs (|17[) and (|18[) . Finally the "running" of the spectral index is given by 

d ti 

n r = — — = Wer/- 24e 2 - 2<, (24) 
d inj 

where 

c = Mpi wm ' (25) 

is another slow-roll parameter. At the horizon crossing epoch, the derivative of H with respect to j leads to 

mm 

d 4 = (26) 

Therefore, when changing j into j(l — aj), we get an approximative expression for H as a function of the 
modified momentum 
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It is obvious that GUP seems to enhance the Hubble parameter so that H(a = 0)/H(a ^ 0) < 1. 

One of the main consequences of inflation is the generation of primordial cosmological perturbations [6!| and 
the production of long wavelength gravitational waves (tensor perturbations). Therefore, the tensorial density 
perturbations (gravitational waves) produced during the inflation era seem to serve as an important tool helping 
in discriminating among different types of inflationary models [70j . These perturbations typically give a much 
smaller contribution to the cosmic microwave background (CMB) radiation anisotropy than the inflationary 
adiabatic scalar perturbations [7lj . 

The tensorial and scalar density fluctuations are given as 
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respectively. Then, the ratio of tensor-to-scalar fluctuations, pt/p s , [66|, IZfJ, ll2| reads 
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In the standard case, this ratio is assumed to linearly depend on the inflation slow-roll parameters [6£ 
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It is apparent that Eq. (|27[) gives an estimation for H in terms of the wave number j, 
start with the equation of motion for the scalar field, the Klein-Gordon equation, 

4> + 3H<j) + dfVty) = 0. 



To estimate 6, we 



(32) 



The 0-term has the same role as that of the friction term in classical mechanics. In order to get inflation from 
a scalar field, we assume that Eq. (|32|) is valid for a very flat potential leading to neglecting its acceleration 
i.e., neglecting the first term. Some inflationary models introduce the slow-roll parameter r]H = —<t>/H(j> = 
—H/2HH. Therefore, the requirement to neglect <j> is equivalent to guarantee that r\H << 1- 
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where the potential itself is model-dependent, for example, V{<f>) = M p i exp [— H a p 0] [73|]. In our model, 
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Then, the tensor-to-scalar fluctuations ratio reads 
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Fig. [T] gives the tensorial density fluctuations p t in dependence on the wave number j (left panel) and on the 
slow-roll parameter e (right panel). In both graphs, a is kept constant, a = 1CP 2 GcV -1 i.e., the upper bound 
is utilized. Also, it is assumed the potential is nearly of the order of the reduced mass M p i i.e., y/2V/M p i ~ 1. 
It is obvious that pt diverges to negative value at low j. Increasing j brings p t to positive values. After reaching 
a maximum value, it decreases almost exponentially and simultaneously oscillates around the abscissa. The 
amplitude of oscillation drastically decreases with increasing j. The right panel shows that j?t(e) oscillates 
around the abscissa. Here, the amplitude of the oscillation raises with increasing e. The oscillation can be 
detected essentially in the CMB spectrum quantizing the primordial residuals of the quantum gravity effects. 



' — ' 1 1 1 1 1 -0.6 1 1 1 1 1 1 

02468 10 02468 10 



Fig. 1: The tensorial density fluctuations p t is given in dependence on the wave number j (left panel) and on slow-roll 
parameter e (right panel), where a is kept constant, a = 10 -2 GeV -1 (lower bound) and it is assumed the y/2V/M p i 
remains constant, ~ 1. 
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Fig. 2: The scalar density fluctuations p s is given in dependence on j (left panel) and on slow-roll parameter e (right 
panel), a and V2V/M p i have the same values as in Fig. [T] 



Fig. [5] refers to nearly the same behavior as that of the dependence of scalar density fluctuations p s on the 
wave number j and e. It is apparent that p s diverges to negative value at low j. Increasing j brings p s to 
positive values. But after reaching a maximum value, it decreases almost exponentially Nevertheless its values 
keep their positive sign. The oscillation of p s (e) is also observed. Here, p s (e) behaves almost similar to Pt{k). 
After reaching a maximum value, it almost exponentially decreases and simultaneously oscillates around the 
abscissa. The amplitude of oscillation drastically decreases with increasing e. 

Fig. Ogives the ratio pt/p s in dependence on e in two cases. The first case, the "standard" one, is given by 
solid curve. The second case, the "modified" case, is given by dashed curve. The latter is characterized by finite 
a, while in the earlier case, a vanishes. Compared to the "standard" case, there is a considerable increase in 
the values of pt/ps with raising e. For the "modifiicd" case i.e., upper bound of a = 1CP 2 GeV -1 , the best fit 
results in an exponential function 
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where /i = 0.875 ± 0.023 and v = 1.217 ± 0.014. For the "standard" case, the results can be fitted by 

(37) 
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which agrees very well with Eq. ([3~T]) . The difference between Eqs. ([55)1 and (f3"?| is coming from the factor in 
the denominator reflecting the correction due to the GUP approach. 

As discussed above, the CMB results and other astrophysical observations strongly make constrains on the 
standard cosmological parameters such as the Hubble parameter H , baryon density rib and even age of the 
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Fig. 3: The dependence of the ratio pt/Ps on the slow-roll parameter e is given in "standard" and "modified" cases. The 
GUP parameter a (in "modified" case) and y/2V/M p i have the same values as in Fig. [TJ The horizontal dashed line 
represents constant ratio pt/ps- 



universe [74ll75| . It turns to be necessary to have constrains on the power spectrum of the primordial fluctuations 

f76j ). This is achievable through the spectral index. From Eq. (f22|) . the scalar spectral index at \[2VjM = 1 
reads 
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The "running" of the spectral index n s is defined by Eq. (|24[) . The results of n r = dn s /d In j arc depicted 
in the right panel of Fig. 21 Early analysis of the Wilkinson Microwave Anisotropy Probe (WMAP) data 
(77l [78j indicates that ri r = —0.03 ± 0.018. As noticed in [78|, such analysis may require modification, as their 
statistical significance seems to be questionable. On the other hand, it is indicated that the spectral index 
quantity n s — 1 seems to run from positive values on long length scales to negative values on short length scales. 
This is also noticed in left panel of Fig. HI where we draw n s vs. u>. Such a coincident observation can be 
seen as an obvious evidence that our model agrees well with the WMAP-data. Recent WMAP analysis shows 
that n s — 0.97 ± 0.017 [llj]. The importance of such agreement would be the firm prediction of inflationary 
cosmology through the consistency relation between scalar and tensor spectra. The physics at the Planck's 
scale is conjectured to modify the consistency relation considerably It also leads to the running of the spectral 
index. For modes that are larger than the current horizon, the tensor spectral index is positive (80j . 
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Fig. 4: Left panel: the spectral index n s is given in dependence on e, where j and a are kept constant (equal 1). 
The "running" of n s is shown in the right panel. The solid curves represent the results from the modified momentum 
j — > j(l — aj) i.e., applying the GUP approach. The dashed curves represent the standard case (unchanged momentum) 
i.e., a — 0. 



IV. THERMODYNAMICS OF QUANTUM STATES IN THE INFLATION ERA 

Based on the general uncertainty principle with the minimal length, the statistical (informational) entropy of 
the FLRW universe described by time-dependent metric is calculated in this section. In section HH the FLRW 
cosmology is briefly reviewed. The minimal length related to the Plank scale can be related to the area of the 
cosmological apparent horizon. The latter would be in turn be related to a. Interestingly, such a relation is 
conjectured to characterize the black holes, where its entropy is proportional to the area of its horizon (l9l. [20j. 
Instead of the event horizon of a black hole, the cosmological boundaries can be chosen to be identical with the 
apparent cosmological horizon. It is assumed that the universe will be in locally thermodynamic equilibrium 
state. Similar to the black hole, the degrees of freedom of a field can be dominant near horizon. 

When assuming that a quantum gas of scalar particles is confined within a thin layer near the apparent horizon 
of the FLRW universe satisfying the boundary condition, the number of quantum states can be calculated using 
the GUP approach. In calculating this, we take into consideration a potential change in the phase space |8l| . 

n(a> ) =*[ dRdedc^ d ?* d P° d P* =±fdR[ ^—^dP R , (39) 
v ' (2tt) 3 J y (1-aP) 4 2tt 2 J J (1-tvP) 4 ' v ' 

where w and Pr are the energy and momentum of the scalar field, respectively. R represents the spacial 
dimension of the layer of interest, where the number of quantum state is to be estimated. As given above, 
a locally equilibrium system is assumed in which the temperature of thermal radiation is slowly varying near 
the horizon, so that the temperature is approximately proportional to the apparent horizon, T oc 1/R. Using 
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natural units c = K = ks = 1, then Eq. (|39|) leads to 
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where ,u is the mass of the field of interest, A = 1 + a 2 [fi 2 + (C - oj)uj] , B = 3H 3 R 3 a 6 uj 3 , C = HPR, 
D = \J —fi 2 + oj(—C + uj). F = fi 2 — uj 2 and / = 1 — R 2 / R 2 A is a function of the comoving time t and di- 
mension of the layer R with Ra = ^/(H 2 + k/a 2 ) 1 / 2 . For D = 1 and because of the small value of a, all high 
orders of a are disregarded. Based on these assumptions the number density of the quantum states is given as 
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where K = y/—f/J 2 + (/ + H 2 R 2 ) uj 2 . The results of n{uS) arc shown in Fig. [5] Left panel shows the results 
at the upper bound of a. The lower bound is given in the middle panel, which would be nearly identical with 
the standard case, where a = 0. Subtracting the results from the lower bound from the results from the upper 
bound seems to result in the affect of the quantum gravity. These results are illustrated in the left panel. 
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Fig. 5: The number density of the quantum states in the inflation era is given in dependence on the energy of the scalar 
field bj. Left panel shows the results at the upper bound of a. The results at the lower bound of a are given in the 
middle panel. The difference between the two bounds is given in the right panel. 
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In Boltzmann limit, the entropy can be directly derived from the free energy is given by — J n{uj) / (exp f3u — 
l)duj. The entropy reads (3 2 f wn(uj) / (Asmh 2 (fiuj/2))duj. The resulting expression is given in Appendix |A"1 The 
results are shown in Fig. [6] In performing these calculations, we set H = 10~ 25 GeV, R = 10 34 GeV~\ / = 1, 
/i = 1 GcV and f3 = 2.4 x 10^ 4 GcV. It has been shown that the thermodynamic first law is fulfilled. When 
the cosmological constant A becomes dominant compared to other forms of matter, the entropy is found to be 
satisfied S = A /4G, where A is the area of the apparent horizon. Again, we distinguish between upper and 
lower bound of a. The lower bound is assumed to be very identical with the standard case. 

Left and middle panels in Fig. [5] show the results at upper and lower bounds of a, respectively. Only real 
values are taken into consideration. The absolute values in the latter case are nearly three orders of magnitude 
larger than that in the earlier case. In the earlier case s(oj) diverges at small to. It shows a kind of saturated 
plateau up to w ~ 2 GeV. This is almost the same as it will be shown in Fig. [7J One of the apparent differences 
is the sign. Here, s(lo) is negative. With increasing uj, it decreases almost exponentially and flips its sign to 
negative one. At lower bound of a, we find that s(lo) behaves almost contrarily. We notice that s(u>) remains 
positive and it decreases almost exponentially. The difference between upper- and lowcr-bound-rcsults is shown 
in the right panel. This difference is assumed to approximately give a qualitative estimation for the effects of 
the quantum gravity on the entropy density. It is obvious that s(uj) remains negative. Increasing u results in 
decrease in the absolute values of the entropy density. It is apparent that negative entropy contradicts the laws 
of thermodynamics. 
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Fig. 6: Entropy density of the quantum states in the inflation era is given in dependence on the energy of the scalar 
field uj. As in Fig. the left and middle panels show the results at the upper and lower bounds of a, respectively. The 
difference between upper- and lower-bound-results is shown in the right panel. 



The dependence of free energy on u> is illustrated in Fig. [7] In the left panel, we show the results at the 
upper bound of a. In doing this, we take into consideration the real values, only. We notice that the free 
energy diverges to negative values at very small values of u. Then, F (u>) makes a plateau up to w ~ 2 GeV. 
With increasing w, the free energy arising from the quantum states switches to positive values. Afterwards, it 
increase, nearly exponentially. The middle panel shows the result at the lower bound of a. We notice that the 
absolute value of a is about three orders of magnitude larger than in the case of upper bound (left panel). Also, 
we notice that lower-bound- values remain negative although they exponentially decay with increasing u>. The 
difference between upper- and lower-bound-results is shown in the right panel. It give a qualitative estimation 
for the effects of the quantum gravity i.e., GUP on the free energy when taking into consideration the quantum 
gravity i.e., applying the GUP approach. 




— w[GeV] 



Fig. 7: The free energy of the quantum states in the inflation era is given in dependence on the energy of the scalar field 
u). As in Figs. [5] and [6] the left and middle panels show the results at the higher and lower bounds of a, respectively. 
The difference between upper- and lower-bound-results is shown in the right panel. 
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V. CONCLUSIONS 

We have studied the effects of GUP on the inflationary dynamics and thermodynamics of the early universe. 
In the first part of this work, we introduce an evaluation for the tensorial and scalar density fluctuations in 
the inflation era. Furthermore, comparing with standard case, the case in which the effects of quantum gravity 
are excluded, we give an estimation for the GUP on all these parameters. The tensorial pt and scalar density 
fluctuations p s arc given in dependence on the wave number j and on the slow-roll parameter e. For a systematic 
comparison, the parameter a is is kept constant, a = 10~ 2 GeV -1 i.e., the upper bound is utilized. Also, it 
is assumed the y2V/M p i ~ 1. We noticed that p t diverges to negative value at low j. Increasing j brings 
Pt to positive values. After reaching a maximum value, it almost exponentially decreases and simultaneously 
oscillates around the abscissa. The amplitude of oscillation drastically decreases with increasing j. Also, pt(e) 
is founf to oscillate around the abscissa. Here, the amplitude of the oscillation raises with increasing e. The 
oscillation can be detected essentially in the CMB spectrum quantizing the primordial residuals of the quantum 
gravity effects. 

The spectral scalar index n s is defined by scalar index. The running of this essential parameter in conjectures 
to shed light on its scaling. The WMAP data indicates that the spectral index quantity n s — 1 seems to run 
from positive values on long length scales to negative values on short length scales. This behavior is confirmed 
in our calculations. The importance of such agreement would be the firm prediction of inflationary cosmology 
through the consistency relation between scalar and tensor spectra. The Planck scale physics is conjectured 
to modify the consistency relation considerably. It also leads to the running of the spectral index, as seen in 
our calculations. For modes that are larger than the current horizon, the tensor spectral index is positive. The 
limitation to the apparent cosmological horizon has been discussed in the present work 

Assuming that a quantum gas of scalar particles is confined within a thin layer near the apparent horizon 
of the FLRW universe which satisfies the boundary condition, the number and entropy densities and the free 
energy arising form the quantum states are calculated using the GUP approach. When taking into consideration 
the quantum gravity i.e., applying the GUP approach, a qualitative estimation for the effects of the quantum 
gravity on all these thermodynamic quantities is introduced. 
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Appendix A: Entropy and free energy 



At temperature T = 1//3, the entropy can be deduced from the number of quantum states, Eq. (|4ip. |8_ 
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sH = p 2 J 



-dui 



4sinh 2 (/3o;/2)' 

18g3 1 fi 2 Q 6 { TJs (-12atanh(a) (/ (5 + a 2 ( M 2 - 3a; 2 )) + 6tf 2 i? 2 a 2 a; 2 ln(R)) + 

H 2 R 2 a 2 u}\n(R)[-9 + 132a; + 2a (240 + 36a; + a (-6^ 2 a; - 3 (2 + ^ 2 ) a; 3 + 12a (41a; 3 + ^ 2 (37 + 41w)) 
4a 3 (a; 4 (51 + 32a; ) - 2^ 2 a; 3 (96 + uj 2 ) + fi 4 (-51 - 96a; + 6a; 3 )))) + 
48a 3 (37 + 34aV) uj 2 ln(w)] ) + 

T 2 (u7 + a ( 60 + 2a ( U + + « 3 (-! + 164a) M 4 + 52a V 6 )) - 

i6/a 2 (6# 4 .R 4 a 3 (11 + 9aV) + f (11 + a (86 + a(-l + 164a)^ 2 + 78a 4 /i 4 ))) - 

3a 4 (l8J/ 6 i? 6 a 3 + / 3 (-1 + 164a + 156a 3 /i 2 ) - \2fH A R A a (11 + 3a 2 (-3 + 4^ 2 ))) u - 



6H 4 R 4 a r (9H 2 R 2 - 2/ (-16 + /i 2 )) a; 2 + 2 (26/ 3 - 72/tf 4 i? 4 - 9H & R 6 ) a 7 uj :i + 

(6H 6 R 6 a 5 a; 2 (18 + w(9 + 2w)) + f (-9 + 480a + 888a 3 (^ 2 + 2a; 2 ) + 136a 5 (3ff + 14^ 2 a; 2 - 2a; 4 )) + 
fH 4 R 4 a 3 (24(33 - 4a;)a; + a 2 (^t 2 (27 + 702a; - 270a; 2 - 8a; 4 ) + a; 2 (162 + a;(405 + 98a;)))) + 
f 2 H 2 R 2 (-9 + 480a + 24a 3 (37^ 2 + a;(33 + 70a;)) + a 5 (3a; 2 (18 + (117 - 62a;)a;) + 8^ 4 (51 + a; 2 ) + 
fj 2 (27 + 2a; (351 + 644a; - 4a; 3 ))))) + 12H 4 R 4 a 5 (-3H 2 R 2 a 2 + 8/ (l + 4a 2 M 2 )) ln(w) + 

-SiyffHRa 2 (f (f + H 2 R 2 ) (-3 + 160a) + 3a 3 (296/ 2 + 15H 4 R 4 a 2 + fH 2 R 2 (88 + 9a 2 )) [i 2 + 
10/ (68/ + 5H 2 R 2 ) a 5 //) In [ fei^/fti, - 2A') / 

(fHRa 2 (/ (/ + H 2 R 2 ) (-3 + 160a) + 3a 3 (296/ 2 + 15H 4 R 4 a 2 + fH 2 R 2 (88 + 9a 2 )) ^t 2 + 
10/ (68/ + 5# 2 A> 2 ) a 5 //) uj)] - 
9/# 3 7? 3 a 5 



V/ + H 2 R 2 



(88 (/ + H 2 R 2 ) + 3 (39/ + 28H 2 R 2 ) a 2 fi 2 ) In [2 (ju) + H 2 R 2 to + sj f + H 2 R 2 K^ J I . (Al) 



We notice that the expression contains complex terms. The results at upper and lower bounds of a are given 
in Fig. [5] (shown in the left and middle panels, respectively). The difference between s(uj) at upper and lower 
bound of a is given in the right panel. Only real values are drawn. 
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Also, the free energy can be deduced, directly, from the number of quantum states, Eq. (|4"Tj) . 

„, n f n(ui) 



i 



— [-12atanh(a) (/ (5 + a 2 (^ 2 - 3a; 2 )) + QH 2 R 2 a 2 uj 2 ln(i?)) + 



im 3 R 2 a & l3 \u) 3 

(-9 + 132a; + 2a (240 + 36a; + a (-6/i 2 a; - 3 (2 + ^ 2 ) a; 3 + 12a (41a; 3 + /i 2 (37 + 41cd)) - 
4a 3 (a/(51 + 32a;) - 2/i 2 a; 3 (96 + a; 2 ) + ^ 4 (-51 - 96a; + 6a; 3 )))) + 
48a 3 (37 + 34aV)a; 2 ln(a;))] + 

T 2 (^3 + a ( 60 + 2a ( 11 + S 6 ")^ 2 + " 3 (-l + 164a) M 4 + 52a V 6 )) - 

^6/a 2 (6# 4 .R 4 a 3 (ll + 9aV) + f (ll + a (86 + a(-l + 164a)^ 2 + 78a 4 /i 4 ))) - 

3a 4 (l8# 6 i? 6 a 3 + f (-1 + 164a + 156a 3 fi 2 ) - 12fH 4 R 4 a (ll + 3a 2 (-3 + 4^ 2 ))) w - 
6H 4 R 4 a 7 (9H 2 R 2 - 2/ (-16 + p 2 )) ^ 2 + 
2 (26/ 3 - 72fH 4 R 4 - 9H 6 R 6 ) aV + 

U + W R \^ K ( 6h6rW ( 18 + -0 + 2 -))+ 
f (-9 + 480a + 888a 3 (f? + 2a; 2 ) + 136a 5 (3// + 14^ 2 a; 2 - 2a; 4 )) + 

fH 4 R 4 a s (24(33 - 4a;)a; + a 2 (^ 2 (27 + 702a; - 270a; 2 - 8a; 4 ) + a; 2 (162 + a;(405 + 98a;)))) + 
f 2 H 2 R 2 (-9 + 480a + 24a 3 (37^ 2 + a;(33 + 70a;)) + 

a 5 (3a; 2 (18+ (117 - 62w)w) + 8/i 4 (51 + a; 2 ) + fi 2 (27 + 2a; (351 + 644a; - 4a; 3 ))))) + 
\2H 4 R 4 a (-3H 2 R 2 a 2 + 8/ (l + 4a 2 /i 2 )) m(w) + 

-iy/fHRa 2 [f (/ + # 2 i? 2 ) (-3 + 160a) + 3a 3 (296/ 2 + 15# 4 i? 4 a 2 + fH 2 R 2 (88 + 9a 2 )) /i 2 + 

10/ (68/ + 5H 2 R 2 ) a 5 /i 4 ] In [(2*v7a* - 2A ") / 
(fHRa 2 (/ (/ + # 2 i? 2 ) (-3 + 160a) + 3a 3 (296/ 2 + l5H 4 R 4 a 2 + fH 2 R 2 (88 + 9a 2 )) ^t 2 + 
10/ (68/ + m 2 R 2 ) a 5 //) w)]- 

y==T=j t 88 C + H2jR2 ) + 3 ( 39 ^ + 2Sh2r2 ) Q V ] In [2 (/a; + H 2 R 2 oj + y/ f + H 2 R 2 K 

Again, we notice that the expression contains complex terms. The results at the upper and lower bounds of a 
are given in Fig. [7J left and middle panel, respectively. The difference between s(oj) at upper and lower bound 
of a is given in the right panel. 
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